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Random walkers absorbing on a boundary sample the Harmonic Measure linearly and indepen- 
dently: we discuss how the recurrence times between impacts enable non-linear moments of the 
measure to be estimated. From this we derive a new technique to simulate Dielectric Breakdown 
Model growth which is governed nonlinear ly by the Harmonic Measure. Recurrence times are shown 
to be accurate and effective in probing the multifractal growth measure of diffusion limited aggrega- 
tion. For the Dielectric Breakdown Model our new technique grows large clusters efficiently and we 
are led to significantly revise earlier exponent estimates. Previous results by two conformal mapping 
techniques were less converged than expected, and in particular a recent theoretical suggestion of 
superuniversality is firmly refuted. 
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I. INTRODUCTION 

The steady state distribution of random walking par- 
ticles is governed by Laplace's Equation. As a result 
Witten and Sander's model of Diffusion Limited Aggre- 
gation (DLA) Q, in which a cluster grows by irreversible 
accretion of dilute diffusing material, has been of dou- 
ble interest: it is readily simulated out to huge cluster 
sizes 0, 1^, whilst at the same time the governance of 
its growth by Laplace's equation renders it a landmark 
mathematical challenge to analyse. The connection to 
a Laplacian field also underpins the breadth of appli- 
cation of the DLA model, to problems such as viscous 
fingering in porous media H| and electrodeposition |^ |^ . 
This interest has all been abetted by controversy as to 
whether the distribution of cluster shapes conforms to 
simple fractal scaling (see and references therein), and 
interest in the multifractal scaling of the growth measure 
[111(1^,111. 

Physical analogies and the mathematical connections 
have led to interest in other models where the growth is 
governed non-linearly by a Laplacian field. In particu- 
lar Niemeyer, Pietronero and Wiesmann introduced the 
family of Dielectric Breakdown Models (DBM) |lj| where 

Vn OC |a„(/)|'' , V^(/) = 0, (/"interface ~ (1) 

and T] is (for interest) a positive parameter. Interest in 
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this is further prompted by Ball and Somfai's proposal 
Is","]!] that growth proportional to field with non-trivial 
spatial cutoff maps onto simple DBM but with rj ^ 1. It 
is important that diffusion controlled growth which is 
locally limited by the capillary energy associated with 
surface ramification is a case in point. Computationally 
these non-linear models have been a challenge, as ran- 
dom walkers only directly sample the harmonic measure 
linearly, realising only the 77 = 1 case. 

First in this paper we show that random walkers can 
be exploited to sample non-linear moments of the har- 
monic measure. In this way we obtain results for the 
active portion of the multifractal spectrum of DLA far 
beyond existing results. The key strength of these meth- 
ods is that no explicit solving of the Laplace equation is 
involved. 

Wc then exploit this to establish a new method of grow- 
ing DBM clusters by random walker accretion. This also 
entails adopting the noise reduction strategies lately in- 
troduced in Ref . 0| , and enables us to explore the DBM 
class out to unprecedentedly large clusters. We show that 
in two dimensions this largely resolves how the exponents 
of the DBM model depend on 77: in particular superuni- 
versality of the tip scaling exponent alpha is strongly 
refuted, in favour of a continuous variation of exponents 
which also confirms the hypothesised upper critical value 

dliaiillii ?7c = 4. 

Our walker-DBM results are supported by extensive 
computations using established iterative conformal map- 
ping methods due to Hastings and Levitov (HL) |23| 
and also by direct integration of the Shraiman-Bensimon 
Equations |i21i] exploiting the mappings of Ball and Som- 
fai [1^ IT3 |. For exponents, all three agree within sta- 
tistical errors. Below the level of the errors there is a 
systematic difference between walker-DBM and HL, and 
separate results for the relative penetration depth sug- 
gest that it is the walker-DBM clusters which are more 
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converged to asymptotic behaviour. Unlike HL and our 
Shraiman-Bensimon integrations, the new walker-DBM 
technique is not Umited to two dimensions and so the 
way forward appears open to a fuU exploration of the 
DEM class in three dimensions. 

II. RANDOM WALKER SAMPLING 

We consider first the problem of sampling the harmonic 
measure of an equipotential surface. The harmonic mea- 
sure is given on the surface by ^ = ^ where outside 
the surface the scalar field 4> obeys a Poisson equation 
— = S{x) with S{x) the source density; typically 
we will be interested in cases where the source is con- 
centrated at points, particularly oo. We can sample this 
measure by introducing random walkers at the source 
points, and tracing their Brownian trajectories to the 
point of first encounter with the surface, whereupon any 
given walker is discarded. The points of first encounter 
uniformly sample the measure fi. 

By firing a large enough sample of random walkers and 
collecting frequency counts of their hit distribution, one 
can build up an approximation to the entire harmonic 
measure. This has been successfully used by Somfai et 
al. However such methods are expensive and give 

differing quality estimates across the support. 

Here we focus on recurrence times defined as follows. 
We first divide the support into (many) small partitions 
(hereafter termed "sites") for each of which we aim to 
estimate the corresponding hit probability /i. We then 
fire (independent) walkers sequentially at the surface and 
when each walker hits, the number of walkers fired since 
the previous hit on that site we will call the "age" a, 
and this provides a simple estimate of the hit probabil- 
ity of that site fxi — 1/a. This is a standard way to 
estimate frequency of uncorrclated events, by their re- 
currence time. The probability distribution of the esti- 
mator 111, given the true underlying value /i for that site, 
is given by simple Poisson statistics as 

Pi{a\fi) = /xe-^°, 

assuming that /i 1 so that we can approximate age 
as a continuous variable. We can generate more reliable 
estimates by using the age ak since the fc'th previous hit, 
which is distributed according to 

Pk{ak\fi)^fie-^^''{fiak)''-^/{k-l}l, (2) 

in terms of which our estimate is fik = k/ ak- 
in applications considered below we will be interested 
in calculating moments of the measure, Mq = X^A**^! 
where the summation is over the sites over which we 
partitioned the support. This can be interpreted as 
Mq = (^'^~^) where the average is weighted by the mea- 
sure ^ itself which in turn is sampled by where random 
walkers hit. At each impact we can use the site age in 
our estimate for the factors /i'""'^, leading to the estimate 

M,fc = ^(fc,(7)((afe/fc)i-«). 



where the average is now over all random walkers hitting 
the surface, we use fc'th ages, and A{k, q) is a trivial 
numerical front factor. 

Random walkers being cheap, our main concern is 
when this estimate converges. Using the distribution (0) 
one can readily check that Mqk converges to Mq pro- 
vided fc > q — 1, and we use the appropriate prefactor 

A{k,q) = ^ i^k-q)]^^' ■ Thus we can compute moments 
of finite degree simply by using high enough order ages. 
In practice we will see that for problems of interest the 
most significant limitation comes not at high q but rather 
at highly negative q, for which random walkers give an 
inefficient sample of the relevant parts of the measure. 

III. MOMENTS UNDER GROWTH 

The simple ideas above become rather useful when ex- 
tended to compute moments of the measure as a sur- 
face grows. The harmonic measure of Diffusion Lim- 
ited Aggregation provides a well studied (but not en- 
tirely resolved) test case. For a cluster of N added par- 
ticles the conventional multifractal scaling would lead to 
Mq N-'^ii)/D ^ where D is the fractal dimension relat- 
ing radius R to N through N ~ . Summing these 
moments over N with weight A^*~^gives us a partition 
function Z{q,t) which we can estimate (ignoring numer- 
ical prefactors) as 

N 

Z(g,<,W) = ^ai-V-\ 

n=l 

where the sum is now over the particles used to grow 
the cluster and the corresponding ages of the sites where 
they hit. Following the spirit of how Halsey et al. 
generalised the identification of multifractal spectrum, 
we can now identify that T{q)/D separates the values 
t < T{q)/D for which Z{q,t,N) — > cx) as ^ oo from 
the values t > T{q)/D for which Z{q, t, N) 0. 

The above definition does not restrict the behaviour 
of Z{q,t,N) on the locus t = T{q)/D but the simple 
expectation is of a logarithmic divergence with N . Then 
a numerical strategy is to choose t such that Z{q, t, N) — 

N 

Z{q,t,eN) — J2 a^^^n*~^ becomes independent of 

n=€N+l 

N as N becomes large with fixed e. 

To obtain results at g > 2 we have to use higher order 
ages. The first order age is naturally thought of as the 
age of the parent to a given new site and a corresponding 
estimate of 02 is given by the age of its grandparent (the 
parent of its parent) and so on. It is of course a concern 
that the more such generations one looks back for the age, 
the more the cluster geometry will have evolved in the 
meantime: we will see that overall cluster screening cuts 
in below a threshold value of g in a well understood way. 
To alleviate the worry above threshold we have studied 
noise reduced clusters where each particle deposited ad- 
vances the local geometry by only a small fraction A of 
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FIG. 1: The moment spectrum r(g) of the harmonic measure 
obtained over the growth of a single DLA cluster by impos- 
ing the condition Z{q, t, N) - Z{q, t, iV/4) = Z{q, t, N/4) - 
Z{q,t,N/16) to yield t = T{q)/D. The respective curves used 
ages of order k = 1,2,4,8 and can be seen to break away 
around q = k + 1 as expected from our discussion in section 
irn Note also close agreement with the electrostatic scaling 
law [2^ that t(3)/Z) = 1. These results come from a single 
off-lattice DLA cluster of = 10® particles grown with noise 
reduction factor 

one notional particle size, so that we can be confident 
that the local geometry does not evolve significantly on 
looking back of order 1/A generations. 

Figure ^ shows the moment spectrum T{q) resulting 
from imposing the condition Z{q,t,N) — Z{q,t,€N) = 
Z{q, t, eN) — Z{q, t, e^N) on a single off-lattice DLA clus- 
ter. As expected from our discussion in sectional the 
results using ages of order k break away at q = fc -I- 1. 
An important check is given by the measured value 
t{3)/D = 1.01 compared to the prediction of imity from 
the Electrostatic Scaling Law. 

For more definitive results we have studied a sample 
of 10 clusters out to iV = 10'' using A = 0.1. Even al- 
lowing for the noise reduction which means that of order 
1/A = 10 particles were deposited per local unit of ad- 
vance in the growth, it is unprecedented to probe the 
harmonic measure out to such large scales. Because of 
the volume of age data, it was cumbersome to seek sta- 
tionarity of Z{q,t,N) — Z{q,t,eN) and instead we ex- 
tracted t{q) from the simple scaling expectation that 
Z{q,t',N) - Z{q,t\eN) oc iV*'-*(«) where we simply 
chose t' to assure reasonably uniform weighting of the 
data. We then took T{q) = Dt{q) using D = 1.71, where 
any error in this value is not significant to the accuracy 
of t{q). The resulting f{a) spectrum is shown in Fig. [21 
for ages of order A: = 1, 3, 9, 27 and compared to one ear- 
lier report of the /(a) spectrum which is distinguished by 
carrying full error bars ,24] and one more recent spectrum 
[25j| which claimed better convergence than all previous. 

The success of our data lies in the region q > 1 corre- 
sponding to a < 1, where it passes two important tests 
by meeting the (dashed) tangent lines corresponding to 
Makarov's Theorem[2^ at g = 1 and (for the higher k 



as appropriate) Halsey's electrostatic scaling law|23| at 
q — 3, and in these respects it clearly improves on the 
earlier published data. 

A significant limitation of the f{a) spectrum obtained 
is that it substantially undershoots at large a, but this is 
qualitatively consistent with expectations. It is obvious 
that we cannot expect to probe probabilities to hit a site 
which are smaller than 1/N as these are unlikely to be 
sampled by the N walkers used to grow the cluster. Thus 
we certainly cannot probe the spectrum for a > D. 

A more careful argument suggests that in principle we 
should (in the limit of large enough N) be able to probe 
all the way up to a = Z), as follows. Consider a typical 
site X "born" with hit probabihty i?"", a < D. It is 
expected that the probability for growth within distance 
r of this point varies as p{r) w [r/R)"-, and when this 
vicinity has received walkers the local structure up to 
length scale r will have been completely regrown. This 
requires a number of walkers added to the cluster given 
by 5N[r/R)°' « leading to 5N « i?"r^-". The ar- 
gument is consistent for a < D because it shows that for 
length scales larger than r the regrowth threshold will not 
have been met, validating our retention of the hit prob- 
ability p(r) ~ (r/i?)". From the point of view of site 
X significant reorganisation of the cluster happens first 
on the smallest scales, and requires 5N « by which 
point we can expect to have hit site X itself with non- 
zero probability. For a > D the scenario is of screening 
being dominated by distant growth and hitting the site 
is unlikely. A more detailed discussion of the screening 
dynamics is given by Ball and Blunt 

Figure 13 shows how the /(a) spectrum from first order 
ages develops as a function of size up to TV = 10^. As ex- 
pected the spectrum builds up on the RHS as we increase 
N consistent with our expectation that ultimately we can 
probe as far as a = 13. Quantitatively the convergence 
is poor beyond about a ~ 1.5. 



IV. DIELECTRIC BREAKDOWN MODEL 

A. The random walker model 

Niemeyer, Pietronero and Wiesmann [T^ introduced 
the generalisation of diffusion controlled growth in which 
the local advance velocity is set by fi^. In this model 
T] — 1 corresponds to DLA, while growth with higher val- 
ues of T] was proposed to model the evolution of dielectric 
breakdown patterns. More recently it has been argued by 
two of us jl^ .14] that non-trivial values of rj can also be 
induced when mapping between different types of ultra- 
violet cutoff, even when the underlying growth is strictly 
proportional to fiux. The computational challenge of the 
DBM is that random walkers sample the harmonic mea- 
sure proportional to /i, so for 77 ^ 1 an explicit calculation 
of the measure appears to be required. Relaxation meth- 
ods [T9I ] and more recentlyffor two dimensions) con- 
formal mapping techniques jlTj ] have been used, but the 
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FIG. 2: Multifractal spectrum of [the harmonic measure of] DLA obtained using ages of order k = 1,3,9,27. Also shown 
are earlier results from Ball and Spivack |2^ and Jensen et al. [2fl| . Note how much better the present data agrees with the 
tangent lines representing Halsey's electrostatic scaling law (for q = 3) and Makarov's theorem (for q = 1). The inset shows 
how the older data extends to higher a, which our method cannot probe. 



resulting computations are dramatically more expensive 
than for simple DLA simulation. 



We propose that the following random walker growth 
model is equivalent to DBM. We fire random walkers 
at the cluster, whose first encounter linearly samples fi. 
At each encounter we use the fc'th order age of the site 
hit to give the estimate fik = fc/ofc and accordingly ad- 
vance the growth locally by an amount proportional to 
6m = Aa].~^ , which is interpreted as the mass added 
to the cluster. There are two clear constraints, the first 
is that the local advance must converge in the mean: 

(^alr^^ oc fj.^'^, requiring that fc > 77 — 1. The second 

is that we desire that the fc'th order ages should look 
back no further than a particle size, so individual advance 
steps should be bounded by 1/fc. Naively this requires 
A < 1/fc for ?7 > 1 and A < N^~^ for 77 < 1, where N 
is the total number of walkers, but we will discuss be- 
low how we can greatly improve on these constraints by 
adjusting A dynamically. 



B. DBM with optimal growth steps 

We now consider allowing the growth step prefactors A 
to vary as the growth proceeds. It is convenient to inter- 
pret A = SQ as the charge borne by each walker adding 
to the growth, and the growth is governed by the local 
advance rate (dr/dQ),^ — fi'^ where Q is the cumulative 
charge added. It was noted by two of us 0| that if we 
characterise the extremal tips of the growth as having 
/^tip ~ that the growth law at the tips trivially 

leads to Q w Ji^+vaup g^g generalisation of a standard 
relation for the fractal dimension of DLA, and less triv- 
ially it is expected that l + ijaup = T(2 + r]) which we will 
use below to render some expressions less cumbersome. 
Within all this framework we are now free to choose the 
charge increments per walker SQ to vary systematically 
with growth of a cluster (but not biased by where each 
particular walker hits), with the corresponding mass in- 
crements given by Sm = SQa^T^ . 

We focus on the case 77 > 1 for which the concern lies 
with anomalously low ages which would give large growth 
steps. At given /i = R~", the probability that a single 
sampling of the fc'th age is less than can be found from 
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FIG. 3: Size dependence of the f{a) spectrum of DLA, ob- 
tained using first order ages and clusters grown per Fig. H to 
different sizes. In tfie vicinity a < D the curves systematicaUy 
rise with the range of A'" used, in support of our claim that 
the age method ultimately probes all the way to a = D. The 
curves are based on the slope of In Z{q,0, N) vs InN across 
successive factors of vTo in A up to iV = 10^ (topmost curve). 



Eq. (|2Jl and varies as {R^^ak)'^ for /ia^ ^1. If we grow 
the cluster by a significant fraction of its size then the 
number of walkers sampling this distribution is given by 
{Q/5Q)R^^°'''~" where the first factor is the total number 
of walkers added and the second is the probability of 
one walker landing on a site with fi « Therefore 
the number of times ages below are expected to be 
sampled is given by 



Q/6QRf'^°'^-°'{R- 



(3) 



where we have taken the maximally contributing value of 
a. To find the smallest likely value of we then set this 
expression to unity. Maximising the a dependence gives 
max„ « R-^(k+i) g^j^j convenient 

to similarly substitute Q k, _R'^(^+''', leading to expected 
minimum age given by 



un{afc} « (i?-(^-+i)-^(2+„)5g^ 



i/fc 



(4a) 



The above calculation is valid when it yields an increasing 
function of R] otherwise the scaling assumption in Eq. © 
does not hold, and we have 



min{afc} = 1. 



(4b) 



We have tested the predicted minimum age by looking 
at the extremal ages sampled during DLA cluster growth, 
where we have SQ = 1, 77 = 1 and t(3) = D. Figure^ 
shows scatter plots of all the ages as a function of cluster 
growth in terms of N R^ , using ages of different orders 
and comparing with the result above. 

We can now exploit the expected minimum age to set 
an optimal choice of SQ in DBM growth: we substitute 
minjafc} into the restriction that the growth steps be 



TABLE I: Numerical values of the parameters used to 
grow walker-DBM clusters. The local growth is 5m — 
{Ao/k)N^{ak/kY-\ 



V 


k 


Ao 




0.5 


1 


0.8 


-0.45 


1.5 


10 


0.85 


0.138 


2 


10 


0.8 


0.234 


3 


10 


0.85 


0.32 



suitably limited Sm <^ 1, leading to 
SQ^R 



(';-!)- 



(fc + l)-T(2 + ,;) 



(5) 



Given that we are restricted to fc-|- 1 > 77, the exponent of 
R in the expression for SQ increases monotonically with 
k towards asymptote (77 — l)a,nin as fc — > 00. 

In practice we are inhibited from setting k too large 
because fc-dependent prefactors compete with optimising 
the power law for finite i?, but taking fc = 10 puts us 
quite close to the limit. Also, we cannot use Eq. 
explicitly because the values of T(fc + 1) and r(2 + 77) are 
not known a priori for a generic 77. Instead, armed with 
the knowledge that such an asymptotic power law exists, 
we use the empirical formula SQ = Aok^^^N^ for which 
the parameters are chosen such that the constraints are 
met. For the parameter values of Tabic used in the 
simulations below, the condition Sm < 1 always holds 
whilst the distance looked back by the k'th order age 
only exceeds one particle diameter a few times per million 
walkers and never exceeds two particle diameters. 

Our algorithm is most competitive for 77 near 1, in- 
cluding the important region 1 < rj < 2 corresponding to 
surface tension regularization [T3 . Il4| | , where it can gen- 
erate cluster masses inaccessible by other methods. For 
too small or large 77 (particularly ?7 ^ 3) it becomes less 
advantageous. 

In principle for A: — > 00 the number of walkers required 
to grow a cluster is given by 



N 



dQ 
SQ 



Jj)H-Qmin + »;(Qtip-amin) 



The result is familiar in that for 7; 1 it recovers the DLA 
result (where N is also the cluster mass) . There has been 
recent controversy about the precise value of atip — amin 
through which our method is explicitly sensitive to 77 [I^ , 
although certainly this difference is very small at 77 = 1, 
where < ayp — amin < 0.03 ± 0.005, see Fig. 



C. Numerical results 

As a test of the new walker-DBM method, we com- 
pared the relative penetration depth (the growth zone 
width, normalized by the average deposition radius) with 
measurements on clusters grown by the more established 
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FIG. 4: Scatter plot of the ages of sites hit for DLA: rj — 1, SQ = 1. The data points are decimated for large A'' for clarity. The 
solid lines correspond to the expected scaling of minimal ages obtained from Eqns. (^t-b); these theoretical estimates appear to 
be suitably cautious. 



HL method, for rj — 2. The asymptotic behavior of the 
penetration depth can be considered a sensitive measure- 
ment, as for DLA it caused considerable controversy in 
the past. As seen on Fig.EJ both methods converge to the 
same asymptotic value of the relative penetration depth, 
and of the two the walker-DBM converges slightly faster. 
This can be understood in terms of the higher effective 
noise reduction levels. 

To probe the scaling, we measured the tip scaling expo- 
nent atip and the fractal dimension D for several values 
of rj. On Fig. we compare these results with measure- 
ments obtained with different methods: HL and by inte- 
grating the Shraiman-Bensimon equation; the latter de- 
scribed in more detail in the next section. The data agree 
within the stated uncertainties, showing a smooth mono- 
tonic transition between the extremes 77 = and t] ^ rjc, 
while below the level of significance point-by-point one 
can see some differences. One difference is that for HL, 
atip shows very little variation between 1 < 77 <2 (seem- 
ingly agreeing with an earlier prediction |l3l while 
the walker-DBM does not show this feature. Another 
difference is that for high rj the walker-DBM yields lower 
values for both atip and D than HL. Although for walker- 
DBM it was not practical to obtain values for ?/ > 3, its 
extrapolation is more consistent with rjc = 4. 

V. SHRAIMAN-BENSIMON EQUATION 

Shraiman and Bensimon observed ^ that for DLA m 
two dimensions the evolution of the complex potential re- 
duces to a non-local problem in one space dimension (plus 
time). This is based on the cluster boundary z{0) ex- 
pressed as a complex function of the imaginary part 9 of 
the complex potential, corresponding to cumulative har- 
monic measure or more loosely, charge. In Refs. Q| 
we adapted this to the DBM class, arguing that chang- 
ing from a fixed cutoff scale in physical space to a fixed 
cutoff in charge space could be offset by adjustment of 



the DBM parameter from 770 to 7^1 — arjQ . Then in terms 
of 

K 
k=l 

the dynamics for DBM growth along a channel of width 
2tt with periodic boundary conditions reduces to 



fe K-j 

+ ^^{{^ + m)j -'2k)tpk-ji^j+7n1pTn (6) 
j—l m— 

with tpa — I. The fixed bandwidth limit k < K corre- 
sponds qualitatively to imposing a cutoff of fixed minimal 
charge on the growing tips, and inevitably means that 
the most screened regions of the physical growth are not 
tracked. Scaling arguments lead to the expectation that 
(ipmi'm) cc TO~^+''i(^~^/"^ ioT 1 <^ m <^ K , but give no 
indication of how wide a range of m is required to see 
the scaling behaviour and hence determine the tip scal- 
ing exponent a, and we will see that this is a major issue 
below. 

The trilinear form of the RHS is efficiently evaluated 
by Fourier methods costing of order K\nK per timestep 
for the whole system, and in the results presented here 
we used K = 1023 and a simple predictor-corrector 
timestepping scheme. We set the timestep adaptively 
such that for every -0^ at each timestep: either the pre- 
dicted and corrected updates Stpk agreed within 10% 
or else ipk was being updated by less than 20%. As 
predictor-corrector is a second order method, the result- 
ing worst case error is of order 5%: even this sounds 
generous but as it was applied to the worst case of 1023 
variables the typical precision achieved was very much 
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FIG. 5: Study of the highest probabihty growth sites and 
the tips of DLA clusters. We launched probe walkers onto a 
fixed non-growing cluster and recorded the number of probes 
landed on each site. This enables to calculate the growth 
probability of each site, including the tip (furthest site from 
the origin), a) The ratio of the growth probabilities of the 
highest hit rate site and the tip: it is a slowly increasing 
function of the cluster mass M. The best fit curve is a power 
law with a small exponent: pmax/ptip ~ A/'""?"""'"'''^, with 
Otip — amin = 0.03 ± 0.005. However, one cannot entirely ex- 
clude logarithmic corrections (see inset); in that case Otip = 
amin. b) The number of sites with higher hit rate than the tip. 
This scales as a power of the radius R with the fractal dimen- 
sion of tips as exponent: np>p^^_^ ~ 7?-^*'p ~ M-'^"?''^. The ex- 
ponent is larger, clearly distinct from zero: /tip = 0.38±0.03. 
On all plots the two datasets are standard DLA {A = 1) with 
10'' to 10* clusters depending on size, and noise reduced DLA 
{A = 0.1) with 4200 to lO" clusters. The number of probes 
was set such that the tip of each cluster received 2500 hits. 
The solid lines are best fit for the data points under them; 
dotted lines are extensions towards excluded data points. 



higher, and we chose 20% on the basis of obtaining re- 
sults statistically indistinguishable under refinement. 
Figure IHI^a) shows the observed behaviour of 

w2ij<k l^jl ) vs fc with increasing time for a rep- 
resentative value of r/i . On logarithmic scales the expec- 
tation is a straight line with slope 1 — 1/a and a difficulty 
is immediately apparent, in that two slope behaviour is 
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FIG. 6: (Color online) Comparison of the relative penetration 
depth H = v'{r2) - {r-)^ / (r) for the walker DBM and 
HL methods. For each method the penetration depth was 
measured on 30 DBM clusters, r] = 2, with 3000 probes on 
each cluster. Both methods converge to the same asymptotic 
value, with the walker-DBM converging slightly faster. 



clearly observed. It is clear that it is the slope emerging 
from lower k values which predominates at large times, 
although different slope is preserved over a limited up- 
per range of fc. Figure |HJb) shows results at long times 
for various rji from which it is clear that whilst the high 
fc slopes are insensitive to 77, the (we claim) asymptotic 
slopes at lower fc exhibit a systematic variation. 

From the present data the conclusion has to be that 
the slopes in the lower range of fc give our best estimate of 
1 — l/a and it is these results which are compared (and 
agree) with the results by other methods in Fig. [TJa). 
The relatively constant slopes at high fc are what domi- 
nated our earlier numerical conclusions in |l3l Il4| which 
had a factor of 10 less range in fc, and appeared to lend 
support to the claim that the slope, and hence a, might 
be independent of rj. That earlier suggestion is clearly re- 
futed, although we cannot rule out more surprises from 
simulations at decades larger K. 



VI. CONCLUSIONS 

We have shown that recurrence times between random 
walkers provide estimates of the local harmonic measure 
which are highly effective for estimating moments and for 
simulating non-linear growth models. 

For the scaling of multifractal moments this technique 
is limited to the regions of the growth which are suffi- 
ciently active that nearby hits dominate the screening 
of a given site, but in that "active range" our results 
for DLA in two dimensions clearly surpass the quality of 
previous data and offer more conclusive support to the 
theoretical exponent relations of Makarov and Halsey. 
Given that our method is also applicable in higher dimen- 
sions, we believe we have established it as the technique 
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FIG. 7: (Color online) a) The tip scaling exponent otip 
and b) the fractal dimension D as a function of rj. The 
four data sets plotted are the walker-DBM, our HL mea- 
surements, earlier HL measurements of Ref. (l^ . and our 
Shraiman-Bensimon results. Because we have to convert the 
SB results back from measurements at fixed ryi to the corre- 
sponding values of rjo = rji/a, these points have error bars 
on both a and 770- The "walker" data point for = 1 is 
in fact a DLA measurement (with very small uncertainties) 
from Ref. ^3] together with the relation atip = D — 1. The 
dotted lines correspond to the limiting behaviors: dense 2- 
dimensional growth (77 — 0) and non-branching 1-dimensional 
growth (?7 > r]c). 



2 4 6 " 

FIG. 8: (Color online) Direct integration of the Bensimon- 
Shraiman equation Q as described in section [3 leads to 
the prediction that i'cum(fc) = Xij<fc IV'jl^ should vary as 
^"jiCi-i/a) ^Yie displayed graphs should exhibit slope 1 — 
1/a. (Recall that rji = ai].) The upper panel a) shows 
plots from the same simulation at rji = 0.5 for successively 
later times (bottom to top). There is characteristic slope 
at high k which persists over a limited range, but it is the 
lower slope which emerges and eventually dominates the wider 
range down to the smallest range of k, and which we take to 
reveal the true tip scaling exponent a. The lower panel b) 
shows data at large times for 771 = 0.3,0.6,0.9,1.2,1.5,1.8 
(top to bottom at right). Whilst the high k slope is insensi- 
tive to 771 there is a clear variation of the asymptotic slopes 
from which we take a in Fig. [T^a). 



of choice to determine harmonic measure in the active 
range. 

Our walker-based simulations of the Dielectric Break- 
down Model have greatly consolidated knowledge of how 
DBM exponents vary with the nonlinearity parameter rj 
in two dimensions. In particular we can now be much 
more confident of the continuous variation of fractal di- 
mension with ?7, looking at the trend in terms of atip 
sharpens the issue, and we can be fairly confident that 
DBM becomes trivial for 77 > rye = 4 as conjectured the- 
oretically [Sllllillil. 

The clearest results for the scaling of the DBM model 
in two dimensions over the full range of 77 now come, iron- 
ically, from our integration of the Shraiman-Bensimon 
equation. However unravelling the confusing two slope 



behaviour which this approach exhibits for the extrac- 
tion of the exponent atip , would not have been prompted 
without the walker-DBM results. Overall it is the agree- 
ment of both methods with earlier results which establish 
the definitive picture. 

Finally we note that the most crucial role of walker- 
based DBM will be in three dimensions, where the tech- 
niques based on complex analysis have nothing to of- 
fer and relatively little is known about the exponent be- 
haviour. We look forward to exploring this in a subse- 
quent paper. 
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